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Abstract. - The particle and current densities are shown to display damping and undergo deco- 
herence in ideal quantum gases. The damping is read off from the equations of motion reminiscent 
of the Navier-Stokes equations and shows some formal similarity with Landau damping. The de- 
coherence leads to consistent density and current histories with characteristic length and time 
scales given by the ideal gas. 



Introduction. — Realistic, complex systems usually 
display separate length scales. In fact, a shorter length 
scale characterizes the elementary excitations, e.g. the av- 
erage separation of particles. The interactions generates 
another, longer length scale, such as the mean free path. 
These scales are separated by some power of the coupling 
constants and the measurements, carried out by macro- 
scopic devices, uncover the manifestations of the longest 
length scale, c.f. collective phenomena. 

One usually thinks of the shortest, microscopic scale 
as the reflection of the dynamics of some trivial, non- 
interacting particle system. It is shown in this letter that 
contrary to this view the observables used to diagnose 
the system may create strong, order one effective inter- 
action vertices for undistinguishable, noninteracting par- 
ticles. What happens is that the correlations generated 
by the exchange statistics of the these particles produce 
involved correlations for local observables and create the 
illusion of some interactions. 

The case of a quantum ideal gas is consedered in this 
work where the effective dynamics of the particle and cur- 
rent densities is sought in the absence of any interaction. 
It is well known that for non-interacting particles these op- 
erators have one-loop, connected Green functions in arbi- 
trary high order in contrast to the elementary fields whose 
higher order Green functions factorize according to Wick's 
theorem. The non-factorization of the Green functions in- 
dicates entanglement and suggests that the effective dy- 
namics of the particle and the current density as local 
degrees of freedom might be highly nontrivial. The main 
result of this work is that decoherence and damping arise 
in the effective dynamics of the particle and current den- 



sity, inherent in an ideal quantum gas. 

Model and method. — Let us cope first with the is- 
sue of interactions in a free system, described by a coordi- 
nate x and the action S[x]. Suppose that we are interested 
in the dynamics of a collective coordinate y = F(x). It is 
easy to see that the extended action 

S[x] -► S[x, y] = S[x] -K j dt[y(t) - F(x(t))} 2 (1) 

with K being a constant generates the correct dynamics 
of both coordinates. Note that the second term in the 
right hand side contains higher than quadratic powers of 
x if the collective coordinate y is a non-linear function of 
x. The action S[x, y] offers another advantage, it allows 
us to consider the free coordinate x as the environment of 
the collective coordinate y. 

We shall consider in this work a model of non- 
interacting particles correlated only by quantum statis- 
tics and investigate the dynamics of the particle den- 
sity and current density, constituting a conserved 4- vector 
j M = (n,j), viewed as collective coordinates (subsystem) 
related to the particle degrees of freedom (environment) 
in a non-linear way, by means of the open time path 
(OTP) formalism. This is a slightly extended version 
of Schwinger's close time path (CTP) method [1] which 
provides the required general framework for exploring the 
open subsystem dynamics. By opening the closed time 
path at the final time we gain access to the density ma- 
trix. 

The CTP formalism has ahead been used to establish 
decoherence and the consistency of histories [2] and to de- 
rive the transport equations [3]. While these approaches 
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aim at the collective mode dynamics arising from gen- 
uine interactions our goal is the construction of the effec- 
tive dynamics of composite operators, particle density and 
current, in ideal gases. A further simplification is that the 
f PI formalism is used below when the experimentally priv- 
ileged local observables, the particle density and current 
operators are considered as the first few terms contribut- 
ing to the McLaurin series of the two point functions in 
the Fourier space. We can make in this manner a short 
cut and consider the effective dynamics of the local density 
and current operators instead of treating the whole com- 
plexity of the two-point functions in the 2PI formalism. 
The result is that we do not have to follow the traditional, 
phcnomenologically motivated way of obtaining the hy- 
drodynamical equations from the energy-momentum con- 
servation but can approach the problem in a more natural 
manner by considering directly the variational equations 
of the effective action for the density and current. 

The path integral representation of the CTP genera- 
tor functional for the connected Green functions of the 
4-current density j M = ip^C^tp [4] reads 



(2) 

where ip and a are CTP doublets and denotes the 
inverse CTP propagator. We use the condensed notation 
where the functions defined in space-time and operators 
acting on them are considered as vectors and matrices, re- 
spectively. The integration over the fields ip and ip' yields 
the non-polynomial influence functional 



W[a] =i^Trln(G- 1 +aC). 



(3) 



The resulting connected Green functions with an arbitrary 
large number of external legs correspond to the one-loop 
structures build up by a single particle line, mentioned in 
the Introduction. We assume that the fluctuations of the 
current are small and retain the quadratic term 



J gI a — —aGa, 



(4) 



where Jg* = i£hTr[GC^] = {<m ,0), n being the parti- 
cle density in the rest frame. The particle-hole propagator 



G^ )(aV) = ithTr[& r ' ,r Cg&" r 'c£]. 



(5) 



contains the exchange statistics factor £. Exploiting trans- 
lational and rotational invariance one finally ends up with 
the CTP structure 



G 



L iS~ 
-iS~ -L 



1 1 
1 1 



where each block X = L, S is a tensor 
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LX t,q 



TX u,q 



(6) 



(7) 



where L = V <g> V/A, T = 1 - L, 



V? =n s hP [ n k 



k,q 



k,q 



hq 2 
2m 



Tlkq 



UJ - 



hq 2 
2m 



with/ k = /rf 3 fc/(2^) 3 , 
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.T&(k+3) 



k,q 



=F— (k - 

1 m \ 



(8) 



(9) 



2, ^(k+§)^(k+f 

is the CTP generalization of the Lindhart function, 5 ± = 
R+±R-, 



w,q 



ttTi J 5(±uj - w k+q + w k )n k (l + £rik + q)F£^, 

(10) 

The retarded prop- 



Bare dynamics. — The bare action 5 s [J] for the 
current expectation value J — SW[a]/Sa, defined by 



n k = l/(e^ - 0, n s = (2s + 1) [5] 
tor is given by G r = —L + iS~ . 



S B [J] + ^Ja 



reads in the quadratic approximation 

S B [J} = l(J~ J P )G-\J- 4). 



(11) 



(12) 



The OTP formalism is realized by sending the final time 
of the underlying CTP scheme to infinity and considering 
this action which contains no implicit boundary conditions 
at the final time anymore for arbitrary, not necessarily 
closed pairs of trajectories (J + , J~). Parametrizing the 
source as a = a(l ± k)/2 ± a [4] and expressing J by 
its CTP components which couple to a (physical current) 
and a (auxiliary current) the current-source coupling aJ = 
a J + aJ a gives J ={,]+ + J-)/2 + n(J+ - J~)/2, J a = 
J + — J~. The inversion of the last two equations yields 
the form of the original current CTP components J ± = 
J — (k=P 1) J a /2 in the bare action in terms of the physical 
and the auxiliary fields J M and J aM . The calculation of 
G _1 of the bare action is tedious but straightforward and 
one finally obtains S B = KS B + i$S B , 



MS B = - 



oj,k"-w,k 



{s tt 



jT :T* iTa 
w,£vw,kJu;,k 



jqT--T* -Ta 



,k'"u,,k' 



,k 

s T 



(f+z2)[(i« fe) 2 + (5 «- ) 



-f iTa* iTa 



i.Ll k f + {S T - k f 



(13) 



with z = muj/hk glis k where fc gas is the characteristic wave 
vector of the one-particle distribution function n k , i.e. the 
Fermi wave- vector hp, for fcrmions and the thermal wave- 
vector \fmk~BT ''/% for bosons without condensate. 
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Consistency and decoherence. The usual way 
of discovering decoherence is either to approximate the 
system-environment interations as a successive scattering 
process [9, 10] or to eliminate a direct product factor of 
the Hilbert space in some simple model [11, 12]. The for- 
mer approach offers a suggestive view of the dynamical 
origin of decoherence but remains rather qualitative. The 
latter, more formal method orients our interest towards 
the degeneracy as the source of decoherence. But to es- 
tablish a sufficiently degenerate environment one needs a 
large number of degrees of freedom and such a system 
can be handled by means of quantum field theory only 
This is the point where the OTP/CTP schemes arise in 
a natural manner because they can handle the (reduced) 
density matrix by means of Green functions in many-body 
systems. 

QS B , which is equal to the imaginary part of the influ- 
ence functional, controlls the consistency [6-8] of the pairs 
of OTP trajectories of the subsystem, running to the fi- 
nal time tf. Furthermore, the dependence of 5sS B on tf 
characterizes the building up of the decoherence [9-12] 
of the environment states |e ± (t/)), generated from |e) by 
the time evolution in the presence of fixed system trajec- 
tories, (e + (t f)\e~ (t /)) = expiS 3 /%. The reduced density 
matrix of the subsystem (ip + \Tr e p(tf)\ip~) = expiW[J]/h 
is given by Eq. (2). If the decoherence of an environment 
state is robust then it is referred to as pointer state [11]. 
A necessary condition for robustness is the persistence of 
the decoherence for macroscopic times, i.e a sufficiently 
large value of QS B for trajectories with uj — > 0. In the 
present consideration the subsystem is characterized by 
the 4- vector (n^kjt.k) for each fixed wave vector k, and 
^sS B indeed controlls the consistency of pairs of small am- 
plitude oscillatory trajectories. 

To construct a measure for consistency and decoherence 
we choose a time dependent state in which the expectation 

value is a plane WELV6 with. cimplitlldG ?T<tJ,k5 Jcj.k 

and con- 
sider the overlap with states with infinitesimally different 
amplitudes. For a single mode, the bare functional integral 
of Eq. (11) assumes the form of ordinary integrals lead- 
ing to the expectation values (n*n) = —2US tt+ (l + z 2 ), 
(n*n a ) = 2ih(L u + iS u -){l + z 2 ), = ~2HS T + , 

(j T *j Ta ) = 2ih(L T + iS T ~) (the indices uj, k being sup- 
pressed). The inverse of the second moment of J a will 
be chosen as a measure of consistence or decoherence [13]. 
Since (n a *n a ) = (j To *j Ta ) = in the absence of external 
source we use instead the ratios 

(n u,2 h\(n*n)\ |S»+| 

1 ' ~ \(n*n a )\ 2 ) 2(1 + z 2 )[{L u ) 2 + (5"-) 2 ]' 

{d t ? = m T *?)\ = \s^\ 

1 ' W*i Ta )\ 2 2[(l*7» + (st-)2]' i > 

for this purpose. Large D indicates that the quantum 
fluctuations around the trajectory with given, oscillatory 
plane wave expectation value for J are distributed by 
mainly classical probabilities. They decohere and the 




Fig. 1: D u (upper) and D T (lower) as functions of Q and z. 



neighboring trajectories tend to be consistent. One ex- 
pects that classical probabilities are recovered at long dis- 
tances and times. The time and distance scale of the 
restoration of classical probabilities can be estimated by 
the value of 1/ui and l/|k| where these ratios become large. 
In the absence of other dimensional constants than h, the 
particle mass m and the density n of the ideal gas these 
scales are in the order of magnitude of fi-/efc gaa , l/fc ga s 
and fall into the quantum domain. Thus consistency or 
decoherence occur at microscopic scale for ideal gases and 
the quantum-classical crossover remains inaccessible for 
macroscopic devices. 

The ratios given by Eqs. (14) are plotted in Fig. 1. 
They are nonvanishing in two strips, starting from the ori- 
gin of the plane (Q, z) and having slopes ±1/2 and width 
2. The center of the strips is the straight line z = ±Q/2, 
corresponding to the mass-shell condition of the environ- 
ment, Twj = h 2 q 2 /2m. The value of the functions along 
this line is shown in Fig. 2. The most consistent density 
trajectories are around (Q,z) — (2, 1). The classical fea- 
tures might be enhanced in the IR directions by means of 
appropriate interactions. The current trajectories become 
more consistent when we move towards the UV direction. 
This does not imply classical behavior at short distances 
or times, the decoherence and consistency being necessary 
but not sufficient conditions for classical physics. There 
are few particles involved at short space or time scales 
and the quantum fluctuations are important. Furthermore 
they do not allow the recovery of the deterministic classical 
laws. The qualitative difference of D u and D T suggests 
that the scale of the quantum-classical crossover is deter- 
mined by the density rather than the current fluctuations. 
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Fig. 2: D u (upper) and D T , (lower) plotted along the line 
z — Q/2 as functions of Q. 



Eqs. (17) yield 
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The longitudinal component of the current, j L , can be re- 
constructed from these equations by means of the continu- 
ity equation for the current. Note that Eqs. (18) contain 
all higher order derivatives. The only approximation in- 
volved is the linearization of the equation of motion. For 
the sake of simplicity (18) is approximated by 



n k 

Jk 



6 Ok, 

e 3 a 



(19) 



We choose a weak external source 



Equations of motion. — The equation of motion for 
the expectation value of the 4-current expectation value 
J = S$iW{a,a]/da is obtained from the effective action 
T[J] = ?RW{a, a] — a - J whose quadratic part in the fields 
is 



kT[J] = JG n - 1 (G r a + J f 



(15) 



Clearly, the effective action is well defined for k ^ only 
[4] . The corresponding equation of motion 



G r -\J-J gl ) = a 



(16) 



is diagonal in Fourier space. The highly complicated, non- 
polynomial equations simplify considerably in the hydro- 
dynamical limit defined by k = k/k gas <C 1, z -C 1, 



-9°a Q 



T-T 

-9 a J 



a + a z 
a + a z 
b + b. 



m tut 

^^gas k 

m iuj 
7ik gas k 



&kk t 2 

"'gas 



^gas 



+ ^)?\ (17) 



^^-gas k 



L.2 
^gas 



where g° = n s mkF/^ 2 , g T = k 2 F /m 2 . For zero tem- 
perature fermions the numerical constants are ao = 1/2, 
a z = 7r/4, a k k = 1/24, b = 3/2, b z = 97r/4 and 
bkk = 39/128. These equations are reminiscent of the 
linearized Navier-Stokes equation of hydrodynamics. The 
different coefficients appearing for the longitudinal and the 
transverse current introduce different bulk and shear vis- 
cosity effects. 

To see the nature of the dissipative terms appearing in 
Eqs. (17) more clearly we consider a stationary flow in the 
w-direction which depends on the x-coordinate only. For 
the static deviation of the density and the current density 
from the corresponding ground state expectation values 



(2< 2 xt 



)3/2 



(20) 



where the index X stands for or T. It yields a Gaussian 
flow, 

„o„.o _ ^ 



'(2^Lw) 3/2 
2 g T u T 

; 3(2^J^ 



(21) 



with £fi ow = \J l 2 xt + l/3kp. Such a spread of the exter- 
nal perturbation is a collective phenomenon due to the 
presence of kp in its length scale and is reminiscent of a 
diffusive process. 

It is easy to find the dynamical origin of this diffusion, it 
is the well known spread of wave-packets of a free particle. 
In fact, the plane wave states not only serve as a useful 
basis for free particles, they represent states which, one- 
by-one, produce robust expectation values. Any smear- 
ing of the wave function of this state in the momentum 
representation to an ordinary wave-packet which is a lin- 
ear superposition of plane waves without some special fine 
tuning of the relative phases leads to the complete loss of 
localization as the time passes, according to the Riemann- 
Lebesgue lemma. The same spread of the individual plane 
wave components spreads any local observable such as the 
density or current. The external source for the three- 
current, a, can be interpreted as a space-dependent shift 
of the momentum variable, eg. p z — > p z + a T in case of 
the external source (20). The spatial spread of such a lo- 
cal Galilean-boost, shown in the second equation of (21), 
indicates the presence of shear viscosity in the dynamics 
of the three-current. 

One finds similar processes in a collisionless plasma. 
The retarded solution of the collisionless Boltzmann equa- 
tion in the presence of an electric field and the Maxwell 
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equations produces Landau damping, the spread of the 
energy of the electric field over the charges [14]. Though 
there is no collision term in the Boltzmann equation, its 
simultaneous solution with the equation of motion for the 
electric field introduces genuine interactions among the 
charges. The spread, described by Eqs. (21) corresponds 
to more elementary processes than the Landau damping 
because it takes place in a truly noninteracting quantum 
gas. A formal similarity between the two phenomena is 
the dephasing, taking place either in the solution of the 
Schrodinger equation for ideal quantum gases or in the 
Fourier integral representation of the retarded Lienard- 
Wiechert potential for classical plasma [15]. 

One should bear in mind the difference between the 
spread of the wave packet and truly diffusive processes. 
The former is related to the asymmetry of the initial 
and final states in physically motivated bases and arises 
in time reflection symmetric dynamics, governed by the 
Schrodinger equation. The nontrivial, dynamical issue of 
the time arrow, the dynamical breakdown of the time re- 
versal invariance appears in the latter when irreversibility 
is observed. The key element is that the soft, gapless col- 
lective modes have a time scale which is much longer than 
the characteristic time of the observation. The bosoni- 
sation is an exact rewriting of the dynamics of a fermion 
system with conserved particle number and it corresponds 
within the functional formalism to the retaining of all 
Green functions of even order as independent variables 
when the fermion fields are integrated over in the path in- 
tegral expressions. Instead of such an extremely involved 
dynamics our calculation can be considered as an approxi- 
mative bosonisation where one keeps the density and cur- 
rent, the first two terms of McLaurin series of the two- 
point funtion, ignoring the remaining informations in the 
two-point funtion and all higher order Green functions. 
Such a coarse graining retains the long distance, low fre- 
quency part of the dynamics and the resulting effective 
theory applies for scales far in the infrared compared to 
the intrinsic scale of the ideal gas. Such a loss of informa- 
tion related to the strong separation of scales generates 
the dynamical breakdown of the time reversal invariance 
and irreversibility. Had we included enough informations 
in the effective description to extend it to the scale of the 
ideal gas we would have lost the dissipative terms. 

Summary. — It is pointed out in this paper that the 
statistics of undistinguishable and noninteracting particles 
produces complicated, nonpolynomial interaction' vertices 
when the dynamics is diagnosed by composite operators. 
In the physically motivated choice of the density and cur- 
rent operators one finds damping and decoherence in an 
ideal gas. These phenomena arise at the only microscopic 
scale of the ideal gas and differ from the analogous phe- 
nomena of truly interactive systems. Namely, the damping 
leads to irreversibility for processes far infrared compared 
to the intrinsic scale of the ideal gas and the decoherence 
of the current operator is stronger at shorter length and 



time scales. These results provide a better background 
to interpret damping and the classical limit of interacting 
systems. In particular, they clarify the physical content of 
the UV intial conditions when the renormalization group 
method is applied to discover the quantum-classical tran- 
sition. 

* * * 
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